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Abstract 

We consider a dynamic server control problem for two parallel queues with randomly varying connectivity and server 
switchover time between the queues. At each time slot the server decides either to stay with the current queue or switch to 
the other queue based on the current connectivity and the queue length information. The introduction of switchover time is a 
new modeling component of this problem, which makes the problem much more challenging. We develop a novel approach 
to characterize the stability region of the system by using state action frequencies, which are stationary solutions to a Markov 
Decision Process (MDP) formulation of the corresponding saturated system. We characterize the stability region explicitly in terms 
of the connectivity parameters and develop a frame-based dynamic control (FBDC) policy that is shown to be throughput-optimal. 
In fact, the FBDC policy provides a new framework for developing throughput-optimal network control policies using state action 
frequencies. Further, we develop simple Myopic policies that achieve more than 96% of the stability region. Finally, simulation 
results show that the Myopic policies may achieve the full stability region and are more delay efficient than the FBDC policy in 
most cases. 

I. Introduction 

Scheduling a dynamic server over randomly varying wireless channels has been a very popular topic since the seminal works 
by Tassiulas and Ephremides in ll28l and ||29l . These works were generalized to many different settings by several authors in 
the network control field (e.g., JU, HU, US), J20], fl2l], J26), OH, (351). However, the significant effect of server switchover 
time between the queues has been ignored. We consider a parallel queue network with randomly varying connectivity and the 
server switchover time between the queues and study the impact of the switchover time on the system performance. 

Our model consists of two parallel queues whose connectivity is varying in time according to a stochastic process and one 
server receiving data packets from the queues by dynamically adjusting its position as shown in Fig. Q] We consider a slotted 
system where the slot length is equal to a packet transmission time and it takes one slot for the server to switch from one 
queue to the other. A packet is successfully received from queue-i if queue-i is connected, if the server is present at queue-i 
and if it decides to stay at queue-i. Therefore, the server is to dynamically choose to stay with the current queue or switch 
to the other queue based on the connectivity and the queue length information of both queues. To the best of our knowledge, 
this paper is the first to consider random connectivity and switchover times to be simultaneously present in the system. Our 
purpose is to characterize the effect of switchover time on system performance. In particular, we are interested in the impact 
of the switchover time on the maximum throughput region (or the throughput region for simplicity) and to find the optimal 
scheduling policy for the server that stabilizes the system whenever the arrivals are within the throughput region. 

Switchover delay in dynamic server control problems is a widespread phenomenon that can be observed in many practical 
systems. In satellite systems where a mechanically steered antenna is providing service to ground stations, the time to switch 
from one station to another can be around 10ms J4], ll30l . Similarly, the delay for electronic beamforming can be on the 
order of 10/is in wireless radio systems (4), ll30l . Furthermore, in optical communication systems tuning delay for transceivers 
can take significant time (/is-ms) 0, ifTTl . We show in this paper that switchover delay indeed fundamentally changes the 
system characteristics. As compared to the seminal work of Tassiulas and Ephremides in ll29l . the supported rate region shrinks 
considerably, the optimal policies change and novel mathematical approaches might be necessary for systems with nonzero 
switchover delay. 
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Fig. 1: System model. Two queues with randomly varying connectivities (C\ and Cj) and t s = 1 slot switchover time. The 
server is currently connected to queue- 1 and it takes 1 time slot to switch to queue-2. 



2 



Note that the switchover time can be smaller or larger than a packet transmission duration in practical systems. In systems 
where it is less than 1 slot, when the server switches from one queue to another, it usually has to waste the entire slot due to 
synchronization issues. For systems with significant switching times (e.g., vehicular networks with mobile relays), our analysis 
can be used as a starting point while keeping in mind that similar solution techniques will apply. Finally note that some of 
our results, in particular the FBDC policy, and the throughput region characterization in terms of state action frequencies hold 
for more general systems such as many queues with arbitrary switchover times and channel statistics. 

We analytically characterize the throughput region A: The set of all arrival rate pairs (Ai,A2) that the system can stably 
support. We derive necessary and sufficient stability conditions on the arrival rate pairs (Ai,A2) in terms of the connectivity 
parameters for both correlated and uncorrected connectivity processes. For this, we consider the corresponding saturated 
system in which there is always a packet to send in both queues and we formulate a discrete time Markov Decision Process 
(MDP) whose stationary deterministic solutions in terms of state action frequencies provide corner points of the polytope of 
achievable rates, i.e., the throughput region. We develop a frame based dynamic control (FBDC) policy for the original system 
with dynamic arrivals. FBDC policy is based on solving a Linear Programm (LP) corresponding to the MDP solution for the 
saturated system and it is throughput-optimal asymptotically in the frame length. FBDC policy is applicable to many general 
systems and provides a new framework for developing throughput-optimal policies for network control. Namely, for any system 
whose corresponding saturated system is Markovian with finite state space, FBDC policy achieves stability by solving an LP 
to find the stationary MDP solution of the saturated system and applying this solution over a frame in the actual system. We 
also develop simple Myopic policies with throughput guarantees that do not require the solution of an LP and that can be more 
delay efficient than the FBDC policy. We show that the Myopic policy with "one lookahead" achieves at least 90% of the 
throughput region while the Myopic policies with 2 and 3-lookahead achieve more than 94% and 96% of the stability region 
respectively. The mathematical solution technique used for proving the stability of various policies is novel in this paper in 
that it involves utilizing Markov Decision Theory inside the Lyapunov stability arguments. 

Optimal control of queueing systems and communication networks has been a very active research topic over the past two 
decades. In the the seminal paper l28l . Tassiulas and Ephremides characterize the stability region and propose the well-known 
max-weight scheduling algorithm. Later in 11291 . they consider a parallel queueing system with randomly varying connectivity 
and prove the throughput-optimality of the Longest-Connected-Queue scheduling policy. These results are extended to the joint 
power allocation and routing problem in wireless networks in ||20| and 1211 and the optimal scheduling problem for switches 
in 11241 and ll26l . Decentralized and greedy scheduling algorithms with throughput guarantees are studied in |6|, Q, fl4l . 1341 . 
while (8l and lfl8l consider distributed algorithms that achieve throughput-optimality (see J9] for a detailed review). In ifTTI . 
l25l and 11351 the network control problem with delayed channel state information is studied, while flj and |13| investigate 
network control with limited channel sensing. These existing works do not consider the server switchover times. Scheduling in 
optical networks under reconfiguration latency was considered in @, where the transmitters and receivers were assumed to be 
unavailable during the system reconfiguration time. While switchover delay has been studied in polling models in the queueing 
theory community (e.g., 12, lfT2l . Ifl5l . OTl ). random connectivity was not considered since it may not arise in classical polling 
applications. To the best of our knowledge, this paper is the first to simultaneously consider random connectivity and server 
switchover times. 

The main contribution of this report is solving the scheduling problem in parallel queues with randomly varying connectivity 
and server switchover times for the first time. In particular, 

• We establish the stability region of the system using the state action frequencies of the MDP formulation for the 
corresponding saturated system. Furthermore, we characterize the stability region explicitly in terms of the connectivity 
parameters. 

• We develop a frame-based dynamic control (FBDC) policy and show that it is throughput-optimal asymptotically in 
the frame length. The FBDC policy is applicable to more general systems whose corresponding saturated system is 
Markovian with finite state and action spaces, for example, networks with more than two queues, arbitrary switchover 
times and general arrival and Markov modulated channel processes. 

• We develop a simple 1 -Lookahead Myopic policy that achieves at least 90% of the stability region while the Myopic 
policies with 2 and 3-lookahead achieve more than 94% and 96% of the stability region respectively. 

■ We present simulations suggesting that the Myopic policies may be throughput-optimal and are more delay efficient than 
the throughput-optimal FBDC policy in most cases. 
This paper provides a novel framework for solving network control problems via characterizing the stability region in terms 
of state action frequencies and achieving throughput-optimality by utilizing the state action frequencies over frames. 

In the next section we introduce the system model and in Section [Til] we provide a motivating example by analyzing the 
case with uncorrelated channel processes over time. We establish the throughput region in Section [IV] via formulating a MDP 
for the saturated system. We prove the throughput optimality of the FBDC policy in Section [V] and analyze simple Myopic 
policies with large throughput guarantees in Section [VI] We provide simulation results in Section I VIII and conclude in Section 
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Fig. 2: Markov modulated ON/OFF channel process. We have e < 0.5 for positive correlation. 



II. The Model 



Consider two parallel queues with randomly varying connectivity and one server receiving data packets from the queues. 
Time is slotted into unit-length time slots equal to one packet transmission time; t £ {0, 1, 2, ...}. It takes one slot for the server 
to switch from one queue to the other, and m(t) denotes the queue at which the server is present at slot t. Let the stationary 
stochastic process Ai(t), with average arrival rate Aj, denote the number of packets arriving to queue i at time slot t where 
E[Af (t)] < A^ ax , i £ {1, 2}. Let C(t) = (Ci(t) 1 C 2 {t)) be the channel (connectivity) process at time slot t, where d(t) = 
for the OFF state (disconnected) and C\(t) = 1 for the ON state (connected). We assume that the processes Ai(t), A-2(t), C\(t) 
and C%{t) are independent. 

We analyze two different models for the connectivity process C(t): 

Definition 1 (Uncorrelated Channels l\20]l , H22V , H29V ): The process C;(i), i £ {1,2}, is in ON state with probability (w.p.) 
Pi and in OFF state w.p. 1 — Pi at each time slot independently from earlier slots and of the other queue. 

Definition 2 (Correlated Channels gj}, $JJ$, $33$, l[36]l): The process Ci(t), i £ {1, 2}, follows the two-state Markov chain 
(i.e., the symmetric Gilbert-Elliot channel model) with transition probability e as shown in Fig. [2] independently of the other 
queue. 

G-E channel model has been widely accepted in modeling and analysis of wireless systems JT], |[T3l , |33l , ||36l , (37). Note 
that our results and algorithms are applicable to general non-symmetric channel models, but here we present the symmetric 
case for ease of exposition. 

Let Q(t) = (Qi(t),Q2(t)) be the queue lengths at time slot t. We assume that Q(i) and C(i) are known to the server at 
the beginning of each time slot. Let at £ {0, 1} denote the action taken at slot t, where at = 1 if the server stays with the 
current queue and at = if it switches to the other queue. One packet is successfully received from queue i at time slot t, if 
m(t) = i, at = 1 and Ci(t) = 1. 

Definition 3 (Strong Stability): A queue is called strongly stable if : 



In addition, the system is called strongly stable (or stable for simplicity) if both queues are stable. 

Definition 4 (Stability Region): The stability region A is the set of all arrival rate vectors (Ai, A2) such that there exists a 
control algorithm that stabilizes both queues in the system. 



The (5-stripped stability region is defined for some S > as A s = |(Ai, A2 ) | ( Ai +5, A2 + 5) £ A|. A policy is said to achieve 
7-fraction of A, if it stabilizes the system for all input rates inside 7A. A throughput-optimal policy achieves 7 = 100% of 
the stability region. 



In this section we show that there is no diversity gain when the channel processes are i.i.d. over time and that channel 
correlation over time is necessary in order to take advantage of the diversity gain and enlarge the throughput region. Specifically, 
we show that when the channel processes are i.i.d. over time, the stability region is reduced considerably with respect to the 
no-switchover time case, and no policy can achieve a stability region larger than that of the simple Exhaustive or Gated type 
policies. Gated policy is such that the server serves all the packets that were present at the queue at the time of arrival and 
then switches to the other queue. In the Exhaustive policy, the server does not leave the current queue until it empties. 

Assume the channel processes C\ (t) and C2 (t) are as described in Definition Q] We first derive a necessary condition on 
the stability of the system and then show the sufficiency of this condition by proving that gated policy stabilizes the system 
under this condition. 

Theorem 1: A necessary condition on stability is given by: 



The proof for a more general system is given in Appendix A. Since both queues have memoryless channels, for any received 
packet from queue-i, as soon as the server switches to queue i, the expected time to ON state is l/pi- Namely, the time to 
ON state is a geometric random variable with parameter Hence, the effect of i.i.d. connectivity is such that this geometric 
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Fig. 3: Stability region under uncorrected (i.i.d.) and correlated (Markovian) channels with and without switchover time. 



random variable is essentially the "service time per packet" for queue- i. Note that we call the term Xi/pi + X2/P2 the system 
load, p, since it is the rate with which the work is entering the system in the form of service slots. In a multiuser single-server 
system with or without switchover times, with stationary arrivals whose average arrival rates are Xi,i £ {1,2}, and i.i.d. 
service times independent of the arrivals with average service times l/pi,i G {1,2}, a necessary condition for stability is 
given by the system load, p, less than 1 . To see this, the stability region of the polling system with zero switchover times is an 
upperbound on the stability region of the corresponding system with nonzero switchover times. Finally, a necessary condition 
for the stability of the former system is p = \\/p\ + X2/P2 < 1. (e-g-, l32l ). Next we show that the stability condition in (fl} 
is also sufficient. 
Gated Policy: 

Serve all the packets that are present at a queue upon arrival at the queue. 

Theorem 2: Gated policy together with cyclic order of service for the server stabilizes the system as long as p < 1, 
The proof for a more general system is given in Appendix B. It is based on a Lyapunov stability argument over a cycle 
duration. Namely, we let m be the discrete time index for the mth time the server stops for servicing a queue and let T m be 
the time slot number of this server-queue meeting times Let I(m) be the i.d. of the node that the mobile serves at time T m and 
let S(Qi( m )(T m )) be the service time required to serve Qi( m ){T m ) packets at time T m . Under Gated service the server serves 
all Qi( m )(T m ) messages, therefore, S(Qjr m -j(T m )) is the summation of Qi( m ){T m ) independent geometric random variables 
of parameter Pjr m ). We have the following queue evolution: 

2 2 2 T m+1 -1 

J2Qi(T m +l)=J2Qi( T ™)+Yl J2Mt)-Ql(m){T m ) (2) 
i=l i=l i=l T=T m 

where T m+ i = S(Qi( m )(T m )) + 1 with additional 1 due to switchover delay. We use the following linear Lyapunov function: 

L(Q(T m )) = V^i^. (3) 
i=i Vl 

Given the current queue sizes, this Lyapunov function represents the expected amount of service slots needed to serve the 
packets present in both queues. We define the drift over one cycle as 

A(T m ) = E{L(Q(T m+2 )) — L(Q(T m ))\Q(T m )} . 

Using ( TT8b and (01 one can show that the drift over the cycle is negative if 

^ftPm) 2 
^ Pi 1 - P 

To understand the intuition behind this condition, first note that is the expected cycle time in the system in steady state 
(in general the expected cycle time is the total travel time per cycle divided by 1 — p) l27l . Hence, Pjzr denotes the expected 

increase in system work load over one cycle. Therefore, ( 1251 ) argues that if 2~2i=i ^'jj" 1 ^ a lower bound on the expected 
decrease in system work load over one cycle, is greater than the expected increase in system load over one cycle, then the 
system is stable. Therefore, the throughput region of the system is given by 

A = {(Ax,A 2 )| ^ + ^<l|. (5) 

I Pi P2 > 

For the case of two parallel queues, the throughput region of the system without switchover delay analyzed in 
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given by 

A ns =|(Ai,A 2 )|Ai <pi,A 2 <p 2 , Ai + A 2 <pi +p 2 (l - Pi)}- (6) 

These two regions are displayed in Fig.[3]for the case of j>\= Vi= 0.5. The stability region of the system without switchover 
time shrinks considerably when there is switchover delay. Note that for the case of deterministic channels, (C\(t) = C 2 (i) = 
1 , Vt), the systems with or without switchover times have the same stability region Ai + A 2 < 1 fl Therefore, there is a significant 
throughput loss due to switchover delay when the channel processes are i.i.d. over time. Therefore, it is the combination of 
switchover delay and random connectivity that result in fundamental changes in system behavior. 

Remark 1: Note that the results of this section hold for more general systems; namely, for systems with N queues and 
arbitrary switchover times between the queues (switchover time from some queue-i to queue- j given by a constant > 1 
slots). The stability region in this case is given by Xi/pi + ... + ^n/pn < L 

With Markovian channels, it is clear that one can achieve better throughput region than the i.i.d. channels case if the channels 
are positively correlated over time. This is because we can exploit the channel diversity when the channel states stay the same 
with high probability. In the following, we show that indeed the throughput region approaches the throughput region of no 
switchover time case in in ll29l as the channels become more correlated over time. Note that the throughput region in l29l is 
the same for both i.i.d. and Markovian channels under the condition that probability of ON state for the i.i.d. channels is the 
same as the steady state probability of ON state for the two state Markovian channels. This fact can be derived as a special 
case of the seminal work of Neely in ETI . 

IV. Stability Region - Correlated Channels 

In this and the following sections we analyze the system under correlated channels assumption. Assume the channel processes 
C\ (t) and C 2 (t) are according to Definition Q] We analytically derive an upper bound on the throughput region of the system 
(necessary conditions on Ai and A 2 for stability) via analyzing the corresponding system with saturated queues. As we show 
in Section [V] the necessary conditions derived in this section are also sufficient and hence the region established in this section 
is the throughput region of the system. 

When switchover times are non-zero, channel correlation impacts the stability region considerably. In particular, channel 
correlation can be exploited to improve the throughput of the system. Moreover, as e — > 0, the stability region tends to that 
achieved by the system with no-switchover time and for < e < 0.5 it lies between the stability regions corresponding to the 
two extreme cases e = 0.5 and e — > as shown in Fig. [3] 

We start by analyzing the corresponding system with saturated queues, i.e., both queues are always non-empty. Let A s 
denote the set of all time average expected departure rates that can be obtained from the two queues in the saturated system 
under all possible policies that are possibly history dependent, randomized and non-stationary. We will show that A = A s . 
We prove the necessary stability conditions in the following Lemma and establish sufficiency in the next section. 

Lemma 1: We have 

A C A s . 

Proof: Given a policy tt for the original system specifying the switch and stay actions based possibly on observed channel 
and queue state information, consider the saturated system with the same sample path of channel realizations for t G {0, 1, 2, ...} 
and the same set of actions as policy tt at each timXe slot t G {0, 1, 2, ...}. Let this policy for the saturated system be n'. Let 
Di(t),i G {1,2} be total number departures by time t from queue-i in the original system under policy tt and let £)^(t),i G 
{1, 2} be the corresponding quantity for the saturated system under policy tt'. It is clear that lim t _ i . 00 (_D 1 (t) + £) 2 (i))/i < 1, 
where the same statement also holds for the limit of D'^t), i G {1,2}. Since some of the ON channel states are wasted in the 
original system due to empty queues, we have 

D 1 (i)<D[(t), and, D 2 (t) < D' 2 (t). (7) 

Therefore, the time average expectation of Di(t),i G {1,2} is also less than or equal to the time average expectation of 
D'iif), i G {1, 2}. This completes the proof since (O holds under any policy tt for the original system. ■ 
Now, we derive the region A s by formulating the system dynamics as a Markov Decision Process (MDP). Let s t = 
(m(t), Ci(t), C 2 (i)) G S denote the system state at time t where S is the set of all states. Also, let a t G A = {0, 1} denote 
the action taken at time slot t where A is the set of all actions at each state. Let H(i) = C(r)|'_ denote the full history 
of the channel processes until time t. For a saturated system, a policy is a mapping from H(i) to the set of all probability 
distributions on actions at G {0, 1}. This definition includes randomized policies that choose at randomly at a given state s f . 

'Throughput region of a general N-queue polling system with stationary arrivals of rates 6 1, ...,N, i.i.d. service processes of mean service times 
si,...,sjv and finite travel times Dij between queues i and j is given by p = X^i=l ^* s « < 1 (see e.g., 1271 ). Therefore, in the absence of random 
connectivity, finite travel times do not affect the stability region. To see this, considering the system under the optimal Gated Policy, with arrival rates close 
to the boundary of the stability region, the fraction of times the server spends receiving packets dominates the fraction of time spent on travel. 
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A stationary policy is a policy that depends only on the current state. In each time slot t, the server observes the current state 
s t and chooses an action a t . Then the next state j is realized according to the transition probabilities P(j\s,a), which depend 
on the random channel processes. Now, we define the reward functions as follows: 

r 1 (« t ,o t ) = lif s t =(l,l,l) or s t =(l,l,0), and a t =l (8) 
r 2 (st,a t ) = 1 if s t =(2, 1, 1) or s t ={2, 0, 1), and a t =l, (9) 

and ri(s tl a t ) = r2(st,a t ) = otherwise. That is, a reward is obtained when the server stays at an ON channel. We are 
interested in the set of all possible time average expected departure rates, therefore, given some u\,u 2 > 0, define the system 
reward at time t as r(st, at) — airi(st, at) + a 2 r 2 (st, a*)- Th e average reward of policy it is defined as 

t=l 

Given some ct\ , a 2 > 0, we are interested in the policy that achieves the maximum time average expected reward r* = max T r 7r . 
This optimization problem is a discrete time MDP characterized by the state transition probabilities P(j\s, a) with 8 states and 
2 actions per state. Furthermore, under every policy, the underlying Markov chain that describes the system state evolution has 
a single recurrent class plus possibly a set of transient states. Note that we eliminate the policy that switches in all 8 states 
and achieves total average rate. Therefore this MDP belongs to the class of Unichain MDPs 11231 . For Unichain MDPs with 
finite state and action spaces, we can define the state-action poly tope, X, as the set of 16-dimensional vectors x that satisfy 
the balance equations 

x(s; 1) + x(s; 0) =V V P(s\s', a)x(s'; a), V s e S, (10) 



ae{0,l} 



the normalization condition 



and the nonnegativity constraints 



J2x(s;l)+x(s;0) = l, (11) 



x(s; a) > 0, for s <E S, a e A. (12) 



Note that x(s; 1) can be interpreted as the stationary probability that action stay is taken at state s. More precisely, a point 
x £ X corresponds to a randomized stationary policy that takes action a € {0, 1} at state s w.p. 

x ( s~ a i 

Ptaction a at state s) = — — ■ — -, a e A, s G S x , (13) 

x(s] 1) + x(s] 0) 

where S x is the set of recurrent states given by S x = {s E S : x(s; 1) + x(s;0) > 0}, and actions are arbitrary for transient 
states s e S/S x 11231 Theorem 8.8.6]. Furthermore, every policy has a unique limiting average state action frequency in X 
regardless of the initial state distribution l23l Theorem 8.9.3]. Therefore, given any policy, there exists a stationary randomized 
policy with the same limiting state action frequencies 11231 . These facts imply that when searching for the optimal policies, 
one can restrict attention to stationary randomized policies as in ( TT3l for x € X. 

The following linear transformation of the state-action polytope X defines the reward polytope lfl6l : {(ji,r2)|ri = x.ri,r2 = 
x.r2, x € X}, where (.) denotes the vector inner product and ri and r 2 are the 16-dimensional reward functions defined in ([8]) 
and (0. This polytope is the set of all time average expected departure rate pairs that can be obtained in the saturated system, 
i.e., it is the rate region A s . An explicit way of deriving A s is given in Algorithm [TJ 

Algorithm 1 Stability Region Characterization 
1: Given ax, a 2 > 0, solve the following Linear Program 

max a\r\ + a 2 r 2 

subject to x e X. (14) 

2: For a given a 2 /a\ ratio, the optimal solution (rj,^) of the LP in (TBI) gives one of the corner points of A s . Find all 
possible corner points and take their convex combination. 



The following lemma is useful for finding the solutions of the above LP for all possible a 2 /ai ratios ll23l Corollary 8.8.7]. 
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Lemma 2: Suppose x is a vertex for the LP in H14[ , then the stationary randomized policy corresponding to x, as defined in 
H13j , is a deterministic policy. Conversely, for any stationary deterministic policy, the stationary distribution of states induced 
by the policy is a vertex for the LP in M4\ . 

The intuition behind this lemma is as follows. For simplicity assume all states are recurrent. Note that the more general 
case can be argued similarly. Now suppose x € X, i.e., x(s, a) > 0, s £ S, a £ A, and x satisfies all the equality constraints 
in ( [Tol l and ( fTTT i out of which only 8 are linearly independent. For a 16 dimensional vector x £ X to be a vertex, we need 
to have at least 16 linearly independent active constraints at x. If x corresponds to a deterministic policy, then either x(s, 1) 
or x(s,0) has to be zero. This gives at least |5| = 8 more linearly independent active constraints at x, satisfying the vertex 
condition. 

Therefore, the corners of the rate polytope A s are given by stationary deterministic policies. There are a total of 2 8 stationary 
deterministic policies since we have 8 states and 2 actions per state. Hence, finding the rate pairs corresponding to the 256 
deterministic policies and taking their convex combination gives A s . Fortunately, we do not have to go through this tedious 
procedure. The fact that at a vertex of ( [Pil l either x(s; 1) or x(s; 0) has to be zero for each s £ S provides a useful guideline 
for analytically solving this LP. The following theorem, proved in Appendix C, is based on this solution to find the corners of 
A s and then applying Algorithm Q] It is one of key results of this paper characterizing the stability region explicitly. 

Theorem 3: The rate region A s is the set of all arrival rates Ai > 0, A2 > that for e < e c = 1 — \/2/2 satisfy 

eA 1 + (l-e) 2 A 2 



(l-e)Ai + (l + e-e 2 )A 2 
Ai + A 2 



(l + e-OAi + (l-e)A 



2 



and for e > e c satisfy 



(l-e) 2 A 1 +eA 2 

Ax + (l-e)(3-2e)A 2 < 
Ai + A 2 < 
(l-e)(3-2e)Ai +A 2 < 
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(1 


-e) 2 
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3 

4 " 
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~ 2 
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4 " 
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~ 2 


< 


3 

4 " 
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~ 2 
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(1 


-e) 2 



(l-e)(3-2e) 
2 

3 e 

4 ~ 2 

(l-e)(3-2e) 



2 

The stability regions for these two ranges of e are displayed in Fig. [4] (a) an d (b). As e — > 0.5, the stability region converges 
to that of the i.i.d. channels with ON probability equal to 0.5. In this regime, knowledge of the current channel state is of 
no value. As e — > the stability region converges to that for the system with no-switchover time in l29l . In this regime, the 
channels are likely to stay the same in several consecutive time slots, therefore, the effect of switching delay is negligible. 

Remark 2: Stability region expressions in terms of the channel parameter e are for two parallel queues, two-state Markovian 
channels and one slot switching time. However, the technique used for characterizing the stability region in terms of the state 
action frequencies is general. For instance, this technique can be used to find the stability region of systems with more than 
two queues, arbitrary switchover times, and more complicated Markovian channel processes. 
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As we stated before, the corner points of the polytope A s correspond to deterministic policies. We observe that the 4 corner 
points of the rate region polytope for e > e c case are also present in the rate region polyopte for e < e c case. One of the two 
additional corners for the e < e c case, namely, the policy for the corner point {r\,T2) = ((1 — e ) 2 /4, (2 — e)/4), corresponds 
to the following deterministic policy: At queue-1: stay only at (Ci, C 2 ) = (1, 0) and switch in all other states and at queue-2: 
switch only at (Ci, C2) = (1,0) and stay at all other states. The critical point is that this policy decides to switch at queue-1 
when the channels are (Ci, C2) = (1,1), which does not provide a vertex for the rate region for e > e c . Analytically, this is 
because (ri,r 2 ) = ((1 — e) 2 /4, (2 — e)/4) is within the convex combination of the 4 vertices of the rate region for e > e c . 
The intuitive reason behind this is that as e increases, the predictions about the future channel states become less reliable. 
Therefore, as e increases, switching at (Ci,C2) = (1,1) at queue-1 for future ON channel states at queue-2 becomes less 
preferable than a successful transmission available at queue-1 in the current slot. 

V. Frame Based Dynamic Control Policy 

We propose a frame-based dynamic control (FBDC) policy inspired by the state action frequencies and prove that it is 
throughput-optimal asymptotically in the frame length. The motivation behind the FBDC policy is that a policy ir* that 
achieves the optimization in ( Tl4b for given weights ot\ and a 2 for the saturated system, should achieve a good performance 
also in the original system when the queue sizes Q\ and Q 2 are used as weights. This is because first, the policy 7r* will 
lead to similar average departure rates in both systems for sufficiently high arrival rates, and second, the usage of queue sizes 
as weights creates self adjusting policies that capture the dynamic changes due to stochastic arrivals. This is similar to the 
structure of the celebrated max- weight scheduling in l28l . Specifically, divide the time into equal-size intervals of T slots 
and let Qi(jT) and Qi{jT) be the queue lengths at the beginning of the jth interval. We find the deterministic policy that 
optimally solves ( TBl i when Qi(jT) and Q 2 (jT) are used as weights and then apply this policy in each time slot of the frame. 
The FBDC policy is described below in details. 



Algorithm 2 Frame Based Dynamic Control (FBDC) Policy 
1: Find the optimal solution to the following Linear Program 

max. { ,. 1: ,. 2} Qi{jT)n + Q 2 {jT)r 2 
subject to (ri,r 2 ) <G A s (15) 

where A s is the rate polytope derived in Section [IV] 
2: The optimal solution (rj , r\ ) in step 1 is a corner point of A s that corresponds to a stationary deterministic policy denoted 
by 7r*. Apply 7r* in each time slot of the frame. 



Theorem 4: The FBDC policy stabilizes the system as long as the arrival rates (Ai, A2) are within the 5-stripped stability 
region Af where 5(T) is a decreasing function ofT. 

The proof is given in Appendix D. It performs a drift analysis using the standard quadratic Lyapunov function. However, it is 
novel in utilizing an MDP framework in Lyapunov drift arguments. The basic idea is that when the optimal policy solving ( fI31 ), 
7r*, is applied over a sufficiently long frame of T slots, the average output rates of both the actual system and the corresponding 
saturated system converge to r*. For the saturated system, the difference between empirical rates and r* is essentially due to 
the convergence of the Markov chain induced by policy ir* to its steady state, which is exponentially fast in T 11611 . Therefore, 
for sufficiently large queue lengths, the difference between the empirical rates in the actual system and r* also decreases with 
T. This ultimately results in a negative Lyapunov drift when A is inside the <5(T)-stripped stability region since from ( fI31 ) we 
have Q(jT).r* > Q(jT).X. 

The parameter S(T), capturing the difference between the stability region of the FBDC policy and A s , is related to the mixing 
time of the system Markov chain and is a decreasing function of T. This establishes that the FBDC policy is asymptotically 
optimal and that A = A s . Moreover, as also suggested by the simulation results in Section [VTll S(T) is negligible even for 
relatively small values of T. 

The FBDC policy is easy to implement since it does not require the solution of the LP for each frame. Instead, we can 
solve the LP for all possible Q 2 (t) /Qi{t) values only once in advance and create a mapping from the Q 2 (t)/Qi(t) values 
to the corner points of the stability region. Then, we can use this mapping to find the corresponding optimal policy ir* at the 
beginning of each frame. Such a mapping depends only on the slopes of the lines in the stability region in Fig. [4] Therefore, 
these mappings are already available and are given in figures [5] and [6] 

Remark 3: FBDC policy is a generic policy applicable to much more general systems. For instance, it provides throughput 
optimality for systems with more number of queues, with swithcover time from some queue-i to queue- j given by a constant 
dij > 1 slots, and with more complicated Markov modulated channel structures. FBDC can also be used to achieve stability 
for classical network control problems such the one with no-switchover times analyzed in |29l . 
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Fig. 5: Mapping from the queue sizes to the capacity region corner points for the FBDC policy for e < e c . 
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Fig. 6: Mapping from the queue sizes to the capacity region corner points for the FBDC policy for e > e c . 



Remark 4: FBDC policy provides a new framework for developing throughput-optimal policies for network control. Namely, 
given any queuing system whose corresponding saturated system is Markovian with a finite state space, throughput optimality 
is easily achieved by solving an LP in order to find the stationary MDP solution of the corresponding saturated system and 
applying this solution over a frame in the actual system. 

Note that the FBDC policy does not require the knowledge of the arrival rates, the channel statistics or the capacity region. 
The mapping in figures [5] and [6] is given in terms of threshold on e since for two queues these values are available. For 
general networks of many queues and arbitrary switchover times, the corresponding table of mappings from the queue sizes to 
stationary deterministic policies can be obtained by solving the LP in ( fT4t using the queue sizes (Qi(jT), Q2{jT) as weights. 

In the next section we consider Myopic policies that do not require the solution of an LP and provide stability for more 
than 96% of the stability region. Simulation results in Section IVIII suggest that the Myopic policies may indeed achieve the 
full stability region while providing better delay performance than the FBDC policy for most arrival rates. 

VI. Myopic Control Policies 

Next, we investigate the performance of simple Myopic policies. We implement these policies in a frame-based fashion 
where the scheduling/switching decisions during a frame of T time slots are based on queue lengths at the beginning of the 
frame and channel predictions for a small number of slots into the future. We refer to a Myopic policy considering k future 
time slots as the fc-Lookahead Myopic policy. In the 1-Lookahead Myopic policy, the server chooses the queue with the larger 
weight where the weight of a queue is the product of the queue length and the expected number of departures in the current 
and the next slot from the queue. The detailed description of the 1-Lookahead Myopic policy is given below. 
Algorithm 3 1-Lookahead Myopic Policy 

1: Assuming that the server is currently with queue 1 and the system is at the jth frame, calculate the following weights in 
each time slot of the current frame; 

Wi(t) = QiO'r)(ci(t) + E[Ci(t + i)|Ci(t)]) 

W 2 {t) = Q 2 (jT)E[C 2 (t + l)\C 2 (t)]. (16) 

2: If W\(t) > W 2 (t) stay with queue one, otherwise, switch to the other queue. A similar rule apply for queue 2. 

Next we establish a lower bound on the stability region of the 1-Lookahead Myopic Policy by comparing its drift over a 
frame to the drift of the FBDC policy. 

Theorem 5: The 1-Lookahead Myopic policy achieves at least ^-fraction of the stability region A s asymptotically in T 
where 7 > 90%. 

The proof is constructive and will be establish in various steps in the following. The basic idea behind the proof is that the 1- 
Lookahead Myopic policy produces a mapping from the set of queue sizes to the stationary deterministic policies corresponding 
to the corners of the stability region. This mapping is similar to that of the FBDC policy, however, the thresholds on the queue 
size ratios Q2/Q1 are determined according to ( TT6b . We first show the derivation of this mapping and then bound the difference 
between the weighted average departure rates of the 1-Lookahead Myopic and the FBDC policies. We refer to the 1-Lookahead 
Myopic policy as the Myopic policy in the following. 

Mapping from queue sizes to actions. Case-1: e < e c For each corner of the throughput region, we will find the range of 
coefficients Qi and Q 2 such that the Myopic policy chooses the deterministic actions corresponding to the given corner. We 
enumerate the corners of the throughput region as 60, •••,&5 where bo is (0,0.5) and 65 is (0.5,0). 



Corner b : 
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Optimal actions are to stay at queue-2 for every channel condition. Therefore, the server chooses queue-2 even when the channel 
state is C\(t),C%(t) = (1,0). Therefore, using (TToT l. for the Myopic policy to take the deterministic actions corresponding to 
b we need 

Qx.(l - e) < Q 2 .(e) => ^ > — . 

Wi e 

This means that if we apply the Myopic policy with coefficients Qi, Q 2 such that Q2/Q1 > (1 — e)/e, then the system output 
rate will be driven towards the corner point bo (both in the saturated system or in the actual system with large enough arrival 
rates). 
Corner b\\ 

The optimal actions for the corner point 61 are as follows: At queue-1, for the channel state 10:stay, for the channel states 11, 
01 and 00: switch. At queue-2, for the channel state 10: switch, for the channel states 11, 01 and 00: stay. The most limiting 
conditions are 11 at queue-1 and 10 at queue-2. Therefore we need, Qi(2 — e) < Q2(l — e) and Q\(l — e) > Q 2 e. Combining 
these we have „ 

2 ~ e < 9l < 1-6 



Note that the condition e < e c = 1 - V2/2 implies that — > 2 



1-e Qi e 

2-e 
1-e ' 



Corner b 2 : 

The optimal actions for the corner point b\ are as follows: At queue-1, for the channel state 10 and ll:stay, for the channel 
states 01 and 00: switch. At queue-2, for the channel states 10: switch, for the channel states 11, 01 and 00: stay. The most 
limiting conditions are 11 at queue-1 and 00. Therefore we need, Qi(2 — e) > ^2(1 — e) and Q\ < Qi. Combining these we 

have _ 

Q2 2-e 

Qi 1-e 

The conditions for the rest of the corners are symmetric and can be found similarly to obtain the mapping in Fig. [7] 
Mapping from queue sizes to actions. Case-2: e > e c 

In this case there are 4 corner points in the throughput region. We enumerate these corners as b , b 2 , b 3 , 6 5 where 6 is (0, 0.5) 
and b 5 is (0.5,0). 
Corner bo'. 

The analysis is the same as the 60 analysis in the previous case and we obtain that for the Myopic policy to take the deterministic 
actions corresponding to bo we need 

Q2 1-e 
Qi > e ■ 

Corner b 2 : 

This is the same corner point as in the previous case corresponding to the same deterministic policy: At queue-1, for the 
channel state 10 and ll:stay, for the channel states 01 and 00: switch. At queue-2, for the channel states 10: switch, for the 
channel states 11, 01 and 00: stay. The most limiting conditions are 10 at queue-2 (since e > e c we have < |5f ) an d 00. 
Therefore we need, Qi(l — e) > Q 2 e and Qi < Q 2 . Combining these we have 

Qi e 

The conditions for the rest of the corners are symmetric and can be found similarly to obtain the mapping in Fig. [8] for e > e c . 
Drift Analysis 

In each frame, FBDC policy drives the system output rate towards the corner point of the throughput region that is the 
solution of the optimization in (fl~5T > (i.e., according to the mappings in figures and |6j. The Myopic policy performs a 
similar operation but according to the different mappings given in figures [7] and [8] In the following we will analyze in which 
Q2(t)/Qi(t) regions the Myopic and the FBDC policies drive the system towards different corner points. We will bound the 
resulting difference between weighted output rates of the two policies where the weights are the queue sizes at the beginning 
of the frames, thereby obtaining a worst case performance for the Myopic policy. Writing the drift expressions for the Myopic 
policy similar to the proof of Theorem |4] we have from d46b that 

<(B + 2)T+Y / Qi(t)Xi-J2 Q l (t)rf Iv +5^2 Q t (t) (17) 
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Fig. 7: Mapping from the queue sizes to the throughput region corner points for the 1-Lookahead Myopic policy for e < e c 
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Fig. 8: Mapping from the queue sizes to the throughput region corner points for the 1-Lookahead Myopic policy for e > e c . 



where r My is the corner point obtained from one of the Myopic mappings and ^(T) is a decreasing function of T. Now let 

W My = J2Qi(t)rf v + (j2Qi(t)\5 3 (T) 

i ^ i * 

and 



W FBDC = X>(*K + X)Qi(*) 



denote the time average weighted departure rates corresponding to the two policies. Also denote the ratio of the two by 
\J> T = WMy/Wopt- We find lower bounds on the ratio ^ over all queue sizes at the beginning of the current frame, which 
will constitute a lower bound on the stability region of the Myopic policy. First, consider the simpler ratio ty' given by 

r _ ZiQi(t)r* Iv 

We later take the 5 factors into account via the argument that iff' — > ty. The following lemma is proved in Appendix E and 
essentially constitutes a lower bound on the achievable throughput region of the Myopic policy. 

Lemma 3: > 0.9002. 

Now consider the expression for ^> given by 

* = 



EiQi(*K + (EiQi(t) js(T) 

Since S(T) and 6s(T) are both decreasing with T, we have that 

* > - 8i(T) 

Choosing T large enough so that S 4 (T) < 0.0002 we have that * > - 0.0002 = 0.9000 and hence W My > 0-9W FB dc- 
Utilizing this in ( fTTb we have the following for the Myopic policy 



^0l<(B+2)T+J2Qi(t)\i-O.9j2Qi(t)r*+S 5 '£Q i (t), 



2T 

% % i 

where 65 = 0.9(5 is a very small and positive number (decreasing with T). Now for (Ai, A2) strictly inside the 0.9 fraction 
of the 85-stripped stability region, there exist a small £ > such that (Ai,A2) + (£,£) = 0.9(ri,r2) — (85,65), for some 
r = (ri,r2) £ A s . Substituting this expression for (Ai, A2) and using J2i Qi(t)( r — r *) < we have, 

2T 



< (b + 2)t- [Y.Qidm 



Therefore, the system is stable for A inside at least the 0.9 fraction of 85-stripped throughput region where £5 (T) is a decreasing 
function of T, 

Remark 5: A similar analysis shows that the 2-Lookahead Myopic Policy achieves at least 94% of A s , while the 3-Lookahead 
Myopic Policy achieves at least 96% of A s . The fc-Lookahead Myopic Policy is the same as before except that the following 
weight functions are used for scheduling decisions: Assuming the server is with queue 1 at time slot t, 

W 1 (t) = Q 1 (jT)(C 1 (t)+J2 k T=1 E{C 1 (t + r)\C 1 (t)}) and W 2 (t) = Q 2 (jT) £^ =1 E{C 2 (i + r)\C 2 (t)}. 
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Fig. 9: The total average queue size for (a) the FBDC policy and (b) the Myopic Policy for T = 10 and e = 0.40. 
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Fig. 10: The total average queue size for (a) the FBDC policy and (b) the Myopic policy for T = 25 and e = 0.25. 
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VII. Numerical Results 

In this section we present simulation results for the FBDC and the 1-Lookahead Myopic policies. We also present numerical 
results that show the stability region for different e values. We performed simulation experiments that present average queue 
occupancy results for the FBDC and the Myopic policies. We first verified the correctness of the simulation model by confirming 
that the FBDC policy achieves the full stability region in the simulation results and then performed experiments for the 1- 
Lookahead Myopic (OLM) policy. In all the reported results, we have (Ai, A2) G A s with 0.01 increments. For each point at 
the boundary of A s , we simulated one point outside the stability region. Furthermore, for each data point, the arrival processes 
were i.i.d., the channel processes were Markovian as in Fig. [2] and the simulation length was 100,000 slots. 

Fig.|9](a) presents the total average queue size, Q avg = Y^l*^ (Qi(t) + Q2(t))/t, under the FBDC policy for e = 0.40 > e c . 
The boundary of the stability region is shown by (red) lines on the two dimensional Ai — A2 plane. We observe that the average 
queue sizes are small for all (Ai, A2) G A s and the big jumps in queue sizes occur for points outside A s . Fig. [9] (b) presents 
the performance of the OLM policy for the same system. The simulation results suggest that there is no appreciable difference 
between the stability regions of the FBDC and the OLM policies. Note that the total average queue size is proportional to the 
average delay in the system through Little's law. For these two figures, the average delay under the OLM policy is less than 
that under the FBDC policy for 86% of all arrival rates considered. For the same system we also simulated a non-frame-based 
Myopic policy that utilizes the queue length information in the current time slot for the weight calculations in ( fToT ). This 
implementation of the OLM policy preserves a similar stability region to Fig. [9] (b) while having delay results at most as much 
as the FBDC policy for 96% of all arrival rates. When the current queue lengths are used in the scheduling decisions, the 
Myopic policy adapts to changes in the system dynamics more quickly, therefore, better delay performance is expected. 

Fig. [TO] (a) shows the total average queue size under (a) the FBDC policy and (b) the OLM policy for e = 0.25 < e c . Again 
this result suggest that the OLM policy is achieving the full stability region. In this case the regular and the non-frame-based 
implementations of the OLM policy outperformed the FBDC policy in terms of delay for 81% and 96% of all arrival rates 
considered respectively. These delay results show that the OLM policy is not only simpler to implement than the FBDC policy, 
but it can also be more delay efficient. 

Figures [TTJ (a) and (b) displays the stability region of the system for e = 0.05 and e = 0.45 respectively. For very small 
channel correlation (e — > 0.5) the stability region tends to that of the i.i.d. channels case, whereas for very large channel 
correlation (e — > 0) the stability region approaches that of the no-switchover time case analyzed in (29). 

Fig. [T2](a) shows the total average queue size under (a) the FBDC policy and (b) the OLM policy for e = 0.10 < e c . Again 
this result suggest that the OLM policy is achieving the full stability region. In this case the regular and the non-frame-based 
implementations of the OLM policy outperformed the FBDC policy in terms of delay for 47% and 91% of all arrival rates 
considered respectively. This suggest that the delay advantage of the Myopic policies are less pronounced for highly correlated 
channels. 
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Fig. 12: The total average queue size for (a) the FBDC policy and (b) the Myopic policy for T = 20 and e = 0.10. 




Fig. 13: The total average queue size for (a) the FBDC policy and (b) the Myopic policy for T = 50 and e = 0.30. 
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Fig. [T3](a) shows the total average queue size under (a) the FBDC policy and (b) the OLM policy for e = 0.30 ~ e c . In 
this case the regular and the non-frame-based implementations of the OLM policy outperformed the FBDC policy in terms 
of delay for 94% and 96% of all arrival rates considered respectively. This suggest that the delay advantage of the Myopic 
policies are more pronounced for less correlated channels. 

VIII. Conclusions 

In this paper, we analyzed the scheduling problem with randomly varying connectivity and server switchover time for the first 
time in literature. We analytically characterized the throughput region of the system using MDP theory, developed a frame based 
dynamic control policy (FBDC) that is throughput-optimal and developed much simpler Myopic Policies achieving 7-fraction of 
the throughput region where 7 can be as high as 94%. For systems with correlated channels, throughput region characterization 
in terms of the state action frequencies of the the saturated system and the throughput-optimality of the FBDC policy hold 
for general systems with many queues, arbitrary switching times and more complicated Markovian channels. Similarly, the 
throughput region as well as the throughput-optimality of the Gated policy for the uncorrected channels case hold for more 
general systems. 

FBDC policy provides a new framework for developing throughput-optimal policies for network control. For any queuing 
system whose corresponding saturated system is finite-state Markovian, FBDC achieves stability based on a novel idea of 
applying state action frequencies that solve an LP for the saturated system. 

In the future, we intend to derive analytical expressions for the throughput regions of more general systems. In particular, 
for systems with non-symmetric Markov channels or multiple-slot switching times, analytical solution of the LP describing the 
throughput region of the system could be possible. We intend to develop throughput optimal Myopic policies for the current 
system and for more general systems. Finally, scheduling and routing in multihop wireless networks with dynamic channels 
and switchover times is an interesting and challenging future direction. 

Appendix A-Proof of Theorem^] 
We prove Theorem Q] for a more general system with A-queues and travel time between queue-i and queue- 7 given by 
slots. We call the term $3 i=1 K/Pi the system load and denote it by p since it is the rate with which the work is entering the 
system in the form of service slots. We prove that a necessary condition for the stability of any policy is p = J^iLi ^i/Pi < 1- 
Proof: Since queues have memoryless channels, for any received packet, as soon as the server switches to queue i, the 
expected time to ON state is l/p%. Namely, the time to ON state is a geometric random variable with parameter p$ and hence 
1/Pi is essentially the "service time per packet" for queue-i. Therefore, the i.i.d. connectivity is essentially a geometric random 
variable representing service time in a classical polling system. In a multiuser single-server system with or without switchover 
times, with stationary arrivals whose average arrival rates are A,,« € {L,2}, and i.i.d. service times independent of arrivals 
with average service times l/pi,i G {1,2}, a necessary condition for stability is given by the system load, p, less than 1. 
To see this, consider the polling system with zero switchover times, stationary arrivals of rate A; and i.i.d. service times of 
mean l/pi- The throughput region of this system is an upperbound on the throughput region of the corresponding system with 
nonzero switchover times (for the same sample path of arrival and channel processes, the system with zero switchover time 
can achieve exactly the same departure process as the system with nonzero switchover times by making the server idle when 
necessary). A necessary condition for the stability of the former system is p = Xi/pi + \n/pn + ■■■ + M/pn < 1, (e-g-, l32l ). 

■ 

Appendix B -Proof of Theorem |2] 

Again we prove the theorem for a more general system with multiple queues and travel time between queue-i and queue-j 
given by Dij slots. We prove that Gated cyclic policy is stable if p = J2 i=1 —■ 

Proof: Let m be the discrete time index for the mth time the mobile stops for servicing a queue. The proof is similar to 
the stability proof in |2). Let T m be the time slot number of this mobile-node meeting times (at time T m+ i the mobile meets 
with the next node in the cycle and at time T m+ N it comes back to the same node). Let I(m) be the i.d. of the node that the 
mobile serves at time T m and let S(Qj( m )(T m )) be the service time required to serve Qi( m )(T m ) packets at time T m . Since 
we have cyclic service, we specify one particular order of service and simplify the notation for traveling times from node i to 
node j, as Di denoting the time required to move from node i to the next node in the cycle. Also let D = J2iLi Di be 
the total travel time in one cycle. 

Since we have gated service, we obtain the following queue evolution: 

N N N T m +S(Q I{m) (T m ))+D I{m) -l 

^Q*(7 m+1 )=^Q ? (T rn )+]r( Mt)) -Qi(m)(T m ). (18) 

i=l i=l i=l t=T m 
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Consider the following Lyapunov function: 



L(Q(T m )) = jr9l^ 



(19) 



The intuition behind this choice of Lyapunov function is that, given the current queue sizes, it is the expected amount of 
service time needed to serve what is currently in all the queues. From dl~8b we obtain, 



Qi{T m+ \) _ ^-v Qi(T m ) 



N T m +S(Q I(m) (T m ))+D Hm) -l 



Pi 



Pi 



E 



i—1 i—1 1 L i—1 

Taking expectations conditional on Q(T m ) we obtain, 



E 

t=T m 



Ai(t)\ Ql(m)(T m ) 



E 



r \ ^ Qi(T m +i) 
p 



Q(T„ 



Qi{T m ) _|_ Aj ^Qj(m)(^m) 



Pt 



Pt 



P/(r 



7(m; 



P/(r 



;j(m 



(^m) 



Pl(* 



(20) 



(21) 



where we used the independence of the arrival and the channel processes conditional on the current queue sizes. 5 , (Q/( m )(T' m )) 
is a random variable that depends on arrivals before T m but not on arrivals after T m as the arrival processes are i.i.d. over 



time. Therefore, E[S(Qi(rn)(T m ))\Qi(m)(Trn)] is nothing but 



Pl( m ) 



Simplifying we obtain 



E[£ Qi{Tm+l) \ci(T m )\ =E 



Qi(T m ) Ql(rn){T m ) 
H P U I{m') ~ I 1 - P)- 



Pl(r 



(22) 



Now we write a similar expression for time T m+ 2- 



E ry Q l (T m+2 ) 

l ~1 Pi 



El E 



N 



Qi{T m +2) 



Q(T, 



m+1, 



Q(T m ) 



w f Qi(^m+l) . n Ql{m+l)(T m+ l) 

= EO h pD I{m+1) (1 - p) Q(T m ) 

I ^ ft Pl(m+1) 

Noting that Q/ (m+ i) (T m+ i) > Qi( m+ \){T m ) and using <E2J, we have from 



(23) 



N 



}i{T m +2) 



Pi 



Q(T m )} < £ + p(D I{m) + D I(m+1) ) - (1 - P ) + 



i=l J ~ 1 »=1 

Repeating the same argument we obtain a drift condition over one cycle given by 



E 



Qi(T m +N) Qi{T m ) 



N-X 



Q(T m ) <pD-(l-p)J2 



!I(m+j) 



(Tm) 



Pl(m+j) 



Hence, we obtain a negative drift as soon as 



N 

E 



Qi(T m ) D 
> P 



1-p 



(24) 



(25) 



Therefore using the Lyapunov stability (e.g., lETl Theorem 3]), the queue length processes at discrete times indexed by m 
satisfies an TV-step negative Lyapunov drift and therefore they are stable. Now consider an arbitrary time slot t G (T m , T„ l+1 ). 
We have that Q{t) < Q{T m+1 ) since there is guaranteed to be no service between T m and T„ l+1 . Therefore we have 
E{Q(t)} < E{Q(T" m+1 )}. Therefore, the system is stable as long as p < 1. ■ 



Appendix C-Proof of Theorem[3] 



We enumerate the states as follows: 



s = (1,1,1) = 1, s = (l,l,0) = 2, « = (1,0,1) = 3, s = (l,0,0)=4, 
s = (2,l,l) = 5 s = (2,l,0) = 6 s = (2,0,l) = 7 s = (2,0,0) = 8. 



(26) 
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We rewrite the balance equations in ( TBI in more details. 

ir(l;l) + x(l;0) = (1 - e) 2 (x(l; 1) + x(5; 0)) + e(l - e)(x(2; 1) + x(6; 0)) 

+e(l-e)(x(3;l)+x(7;0)) + e 2 (x(4; 1) + .t(8; 0)) (27) 
x(2;l)+x(2;0) = e(l-e)(x(l;l) + x(5;0)) + (1 - e) 2 (x(2; 1) + 1(6; 0)) 

+e 2 (a;(3; 1) + x(7; 0)) + e(l - e)(x(4; 1) + x(8; 0)) (28) 

x(5;l)+x(5;0) = (1 - e) 2 (x(5; 1) + x(l; 0)) + e(l - e)(x(6; 1) + x(2; 0)) 

+e(l-e)(x(7;l) + x(3;0)) + e 2 (x(8; 1) + x(4; 0)) (29) 
x(7;l)+x(7;0) = e(l-e)(x(5;l)+x(l;0)) + e 2 (x(6; 1) + x(2; 0)) 

+(l-e) 2 (a;(7;l)+a-(3;0)) + e(l - e)(x(8; 1) + x(4; 0)) (30) 



The following equations hold for each channel state pair (Ci, C2). 

1) + x(l; 0) + x(5; 1) + x{5; 0) = 1/4 (31) 

a:(2; 1) + x{2; 0) + x(6; 1) + x(6; 0) = 1/4 (32) 

x{3; 1) + s(3; 0) + x(7; 1) + x(7; 0) = 1/4 (33) 

x(4; 1) + a:(4; 0) + x(8; 1) + 1(8; 0) = 1/4 (34) 



Let u\ = (x(l; 1) + x(2; 1) and u 2 = (x(5; 1) + x(7; 1)). Summing up d27} with d28j and with (O we have 

e Ul = -(x(l;0)+x(2;0)) + e(x(3; 1) + x(4; 1)) + e(x(7; 0) + x(8; 0)) + (1 - e) (x(5; 0) + x(6; 0)) 
eu 2 = -(x(5;0) + x(7;0)) + e(x(6; 1) + x(8; 1)) + e(x(2; 0) +x(4;0)) + (1 - e)(x(l; 0) +x(3;0)) 



Rearranging and using (f3TT>-(f34T> we have 

ui = + e(x(3; 1) + i(4; 1) + s(7; 0) + x(8; 0)) - (2 - e)(x(l; 0) + x(2; 0)) - (1 - e)(x(5; 1) + a?(6; 1)) 

(35) 

ua = ^^ + e(x(2;0)-x(4;l)+x(6;l)-x(8;0)) - (2 - e)(x(5;0) + x(7; 0)) - (1 - e)(x(l; 1) + x(3; 1)) 

(36) 



Using (|29j in (j35j and ([27]> in d36]» we have 

ui = ^ + e(x(3; 1) + x(4; 1) + x(7; 0) + x(8; 0)) - ^— ^ (x(4; 0) + x(8; 1)) - (1 ~ ^ ~ 2e) x(6; 1) 

+ ^(fef (5; 0) ~ 1 ^2-e) x{1] 0) ~ (X(3; 0) + X(7; 1}) ~ ( 2 ~ £ + ^7") 2:(2; 0) (3?) 

u-2 = + e(x(2; 0) + x(6; 1)) - (e + f£-=i>) (s(4; 1) + x(8; 0)) - (1 = 2&) x(3; 1) 

+ 7(2^7) x{1; 0) ~ l e(2-l) x{5 ' 0) ~ ^7" (a:(2; 1} + x(6; °)) - ( 2 - e + %Z7~) 3; t 7 ; °)- < 38 ) 



IS 



Using (|33j and (04) in $37} and ([32]i in (|38} we have 

Ul = (1 ' 4 (2 (3 e) 2£) + ( e + ^Y^r) 1) + *(8; 0)) + ^ (*(3; 1) + «(7; 0)) 

- U + ^*<* 0) - ^.(1, 0) - (2 - e + ^>(2; 0) (39) 



- i ^ M ^^') + ^)*»)-^*0)-( 2 -« + ^)*0, C40, 

Consider the LP objective function ai(x(l; 1) + x(2; 1)) + a2(a;(5; 1) + x(7; 1)), and note that the solution to this LP is 
a stationary deterministic policy for any given ai and «2- This means that, for any state s either x(s; 1) or x(s; 0) has to be 
zero. In order to maximize ai(x(l; 1) + x(2; 1)) + ot2(x(5; 1) + x(7; 1)) we need 

x{7;0) = if — > 



ai " (2-e)2 + (l-e) 



2 ' 



Oil 1 

x(3; 1) = if — > 



(2-e) 2 + (l-e) 2 > — >1 



ai " (l-e)(3-2e)' 

X (5;0)=0 if ^> 1 ^~ £ 2 , 
ai 1 + e — 

x(8;0) =x(4;l) =0 if — > 1, 

x(6;0)=0 if — > (1- e)(3-2e), 

ft -n n - f "2 . 1 + e - e 2 

xfl; 1) = if — > . 

ai 1 — e 

Note that we have 

(2-e) 2 + (l-e) 2 > (l-e)(3-2e)> 1 

+ e - i 
1 - e 

holding for all e G [0,0.5]. Consider the following two cases: 
Case-1: e > e c = 1 - V2/2 

In this case we have (1 — e)(3 — 2e) < (1 + e — e 2 )/(l — e). This means that we have the following optimal policies depending 
on the value of a?,/a.\. 

l<^<(l-e)(3-2e): 

@queue-l : (1, 1,1) : stay, (1, 1, 0) : stay, (1, 0, 1) : switch, (1, 0, 0) : switch. 
@queue-2 : (2,1,1) : stay, (2, 1, 0) : switch, (2, 0, 1) : stay, (2, 0, 0) : stay. 

Substituting the above zero variables into d39l and d40b . it can be seen that this policy achieves the rate pair 

(l-e)(3-2e) 3 - 2e 

r l = 777; s ) r 2 = 



4(2- C ) ' 4(2 -e) 



si>(l- e )(3-2e): 

@queue-2 : (2,1,1) : stay, (2, 1, 0) : stay, (2, 0, 1) : stay, (2, 0, 0) : stay. 

In this case it is optimal to stay at queue-2 for all channel conditions. Therefore the decisions at queueu-1 are arbitrary. Namely, 
it is sufficient that at least one state corresponding to server being at queue- 1 take a switch decision, which is the case for 
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ocz/oii > ((1 — e)(3 — 2e)), since x(3; 1) = if otz/otx > 1/((1 — e)(3 — 2e)). Since the policy always stays at queue-2, it 
achieves the rate pair 

n = 0, r 2 = 0.5. 
Note that the case for a^/cti < 1 is symmetric and can be obtained similarly. 



Case-2: e < e c = 1 - V2/2 

In this case we have (1 — e)(3 — 2e) > (1 + e — e 2 )/(l — e)- This means that before the state x(6; 0) becomes zero, namely 
for (1 + e— e 2 )/(l — e) < oiz/oti < (1 — e)(3 — 2e), having x(l; 1) = is optimal. This means that there is one more corner 
point of the rate region for e < e c . In more details we have the following optimal policies. 

1 < £2. < — e 2 ■ 
— cti — 1 — e 

@queue-l : (1, 1, 1) : stay, (1, 1, 0) : stay, (1, 0, 1) : switch, (1, 0, 0) : switch. 
@queue-2 : (2,1,1) : stay, (2, 1, 0) : switch, (2, 0, 1) : stay, (2, 0, 0) : stay. 

This policy is the same policy as in the previous case and it achieves the rate pair 

_ (l-e)(3-2e) 3 - 2e 

n ~ 4(2 -e) ' T2 ~ 4(2 -e)' 



22 \ 1+e — e 2 ■ 
ai ^ l-£ - 

We have the following deterministic actions. 

@queue-l : (1,1,1) : switch, (1,1,0): 7, (1, 0, 1) : switch, (1, 0, 0) : switch. 
@queue-2 : (2,1,1) : stay, (2,1,0):?, (2, 0, 1) : stay, (2, 0, 0) : stay. 

In order to find the final threshold on o^/ai, we substitute the above deterministic decisions in (l28l l. ( |29| i and < f30b . Utilizing 
also (EB, OS, (O and (|34> we obtain 

x(2; 1) = - (1 - 6)^(6; 1) (41) 

s(5; 1) + x(7; 1) = ^^ + ex(6;l) (42) 

(43) 

The previous threshold on a^/ai for a;(6;0) to be zero, i.e., (1 — e)(3 — 2e), is valid for the case where x(l;0) = 0. Other 
decisions staying the same, when x(l;0) is positive and x(l; 1) = 0, r 2 increases and n decreases. Therefore the threshold 
on a-i/ax for x(6; 0) to be zero changes, in particular it becomes c^/cti > (1 — e) 2 /e. This gives the following two regions: 

l+e-e 2 < a < (1-e) 2 . 
1 — e — ai — e 



(1,0,1) : switch, (1, 0, 0) : switch. 
(2,0,1) : stay, (2, 0, 0) : stay. 

2 - e 
" 4 ' 



The optimal policy is 

@queue-l : (1, 1, 1) : switch, (1, 1, 0) : stay, 
@queue-2 : (2, 1, 1) : stay, (2, 1, 0) : switch, 

From d4~TT i and d42l it is easy to see that this policy achieves 

(1-f) 2 

»'l = ; ' ? '2 



> 



(1- 



The optimal policy is 

@queue-2 : (2,1,1) : stay, (2, 1, 0) : stay, (2, 0, 1) : stay, (2, 0, 0) : stay. 
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T=0 T=0 



T-l „ T-l 2 T-l T-l 



This policy achives 

ri =0, r 2 = 0.5. 

Similar to Case-1, the case a^/ai < 1 is symmetric and can be solved similarly. 

Thus we have characterized the corner point of the stability region for the two regions of e. Using these corner points, it is 
easy to derive the expressions for the lines connecting these corner points, which are given in Theorem [3] 

Appendix D-Proof of Theorem|4] 

Proof: Let Di(t) be 1 if there is a departure from queue-i at time slot t and zero otherwise, we have the following queue 
evolution relation. 

Qi{t + l) = Q i (t) + A i {t)-D i {t). 

Writing similar expressions for time slots t G {t + 2, ...,t + T} and summing all the expressions creates a telescoping series, 
yielding 

T-l T-l 

Qi(t +T) = Qi{t) - D i(t + T ) + J2 Ai (t + t). 
Taking the square of both sides we obtain 

T-l 

Qi(t + Tf < Qi{tf + ( J2 D i(t + t )Y + (E + T )) z - 2Qi ^ ( D & +t ^-J2 A ^ + t )) • < 44 > 

T = T=0 T=0 T=0 

Define the quadratic Lyapunov function 

2 

L(Q(t)) = 

i=l 

and the T-step conditional Lyapunov drift 

A T (t) 4 E {L(Q(t + T)) - L(Q(f))|Q(t)} . 

Summing (|44| | over both queues, taking conditional expectation, using Di(t) < 1 for all time slots t, E{Ai(t) 2 } < A^ ax and 
EiAiitJAifa)} < Vlp^OpEXA^Jp < A^ ax for all and t a we have 

A T (t) < 2BT 2 +2E I 5^Q<(t)5^ [A ( (t + r) - A(* + t)] |Q(t)| 

I i t=0 J 

= 2BT 2 + 2T]r Qi(t)Xt-2 £ Q 1 (t)E j £ +r)|Q(*)| 

i i I r=0 J 

where = 1 + A 2 nax is a constant. 

Let (t) be a reward function such that fj (t) = 1 if at time i the server is at queue-i with ON channel and decides to stay 
at queue-i at time t and Ti(t) = otherwise. Note that Ti(t) is simply the reward function associated with applying policy it* 
to the saturated queue system whose infinite horizon average rate is r* = (r*,^). Let x* be the optimal vector of state action 
frequencies corresponding to it* . Define the time average empirical reward from queue-i in the saturated system, ?T,i(t), and 
that in the actual system, DT,%{t), as 

?T,m = ^ E r *( <+T )> = i 

r=0 r=0 

Also define the corresponding two dimensional vectors Yx(t) and Dy(t). Similarly define time average empirical state action 
frequency vector xt(£). Let a.b denote the inner product for vectors a and b. From the definition of the rewards in terms of 
state action frequencies in (O we can write Pr,i(i) = a-i-xr(£)> rr,2(i) = a2.Xr(t) and r\ = ai.x*, r\ = a2.x*, where ai 
and a2 are appropriate vectors of dimension 16. Now we have that as T increases, xr(£) converges to x* and hence fr(i) 
converges to r* with probability 1 regardless of the initial state of the system. More precisely we have the following lemma 

EB> ED: 

Lemma 4: For every choice of initial state distribution, there exists constants C\ and C2 such that 

P(||xr(*) - x*|| > <5i) < c ie - C2 ^ T , VT > 1, VSi > 0. 
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Therefore using 1 1 (t) — r*|| < ||xt(£) — x* || (||ai| | 2 + | [ 02 1 1 2 ) 2 , we have that there exists constants c 2 and C3 such that 



P(||f T (i) - r* [J > < cie" C35 i T , VT > 1, V<$i > 0, 



under policy n* for any initial state distribution. Now define the following: 

T-l 



W T (t) 4 ^Q i (i)I^A(i + r)=^Q i (^(i). 

i r— i 

T-l 

R T (t) ± E^w^E^+^E^w^w. 



T = 



We rewrite the drift expression as 



Now we bound the last term. 



Consider 



^1<BT +X)Q i (f)A i -E{W T (f)|Q(t)} 
= BT +5^Q i (f)A i -E{i2*(t)|Q(t)} 

+E {#*(*)- W r (t)|Q(t)} 
= BT +^<9i(^-E^K 

+£{ir(t)-w r (*)|Q(*)}- 

E {£*(*)- Wr(t)|Q(*)} = 

= E {R*(t) - W T (t)\Q(t),R*(t) - W T (t) > <J 2 ||Q(t)||} 
.P(R*(t)-W T (t)>S 2 \\Q(t)\\ |Q(t)) 

+ E {#*(*) - Wr(t)|g(t),ii*W - Wr(*) < *a||Q(t)||} 
.P(R*(t)-W T (t)<S 2 \\Q(t)\\ |Q(t)) 
<(^Q(4p(|i?*(t)-Wr(t)|>5 2 ||Q(t)|| |Q(t))+fc||Q(f)||. 



P(|iE*(t)-Wr(*)|>fc||Q(t)|| |Q(t)) 
<p(V(*)-iM*)l >|||Q(t)|| |Q(i) 



P(|Wr(t)-Rr(t)|>-|||Q(t)|| |Q(t) 



<P(||r*-M*)ll>f |Q(*) 



+ p(||D T (t)-f T (t)||>^|Q(t) 
where the last inequality follows from the Schwartz Inequality for inner products given as 

\R%t)-R T (t)\ = m).(r*-Mm<\\Q(m-\\r*-r T (t)\\- 
Using d45l l. there exists constant C4 such that 

E {R* (t) -W T (t)\Q(t)} <{J2 Q(t))c ie - c ^ T +5 2 ||Q(i)|| 

i 

^Q(t|pf||6 r (t)-f r (*)ll>^|Q(*) 



(45) 
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Hence we can write the drift term as 

^ < bt+ e Qi - E q*(*k + (E 



2T 

2 2 I 

Cie -C45 2 2 T +52+p ^ TW _ fr(t) || >| | ( 

Note that ||Z)y(t) — rr(i)ll is because of the lost rewards due to empty queues and it is equal to zero if both of the queues 
have more than T packets at time t. Namely, ||-D T (i) - r T (t)\\ = if Q\{t) > T and Q 2 (t) > T. Therefore, calling 

8 = c\er Ci&,2 - T + 8 2 , we can write 



A T (t) 
2T 



< 



BT+Y^Qi(t)\i^Qi{t)r*+(5+l {Q(t)<T . 1} %2Q(t) 



< BT+^Q i {t)\ i -^Q i {t)r*+8^Q i (t) + 2T. (46) 

i i i 

Now for (Ai, A2) strictly inside the 8-stripped throughput region A*, there exist a small £ > such that (Ai, A2) + (£,,£,) = 
(fi, f 2 ) — (5, 8), for some r = (n, r 2 ) € A s . Therefore we have, 

^ < (b + 2)t + e - r*) - eE &(*)• 

Finally using ^ - Qi(t)(r — r*) < we have 

A T (t) 
2T 



<(B + 2)T-(E0i(*))f- 



Hence the queue sizes have negative drift when YliQi(t) * s outside a bounded set. Therefore the system is stable for A 
within the 8-stripped stability region Af where 8{T) is a decreasing function of T (see e.g., ETl Theorem 3]). Note that 
8 = c\eT Ci& ?- T + ($2 for any 8 2 > 0. Therefore choosing S 2 appropriately (for example, 8 2 = T~°- 5+Sa for some small 8$ > 0), 
we have that S(T) is a decreasing function of T. ■ 

Appendix E-Proof of LemmaO 

Here we prove that > 0.9002 where 

Proof: We divide the proof into separate cases for different e regions. 
1) Weighted Departure-Rate Ratio Analysis, Case 1: e < e c : Considering the mappings in figures [7] and [5] the regions where 
the Myopic policy and the optimal policy "chooses" the same corner point, we have vp' = 1. In the following we analyze the 
ratio in the regions where the two policies chooses different corner points. We term these cases as "discrepant" cases. We will 
use Qi and Q 2 instead of Qi(t) and Q 2 (t) for notational simplicity. Note that we have that |5f > ^(l-e) always holds. 
However equals (1 ~ £)2 at e t = 0.245 for the case of e < e c = 0.293. 

Case 1.1: e < e t -> f=| < 
Discrepant Region 1: ^ 1 ~ e ^ < ^ < 

In this case the Myopic policy chooses the corner point b\ whereas the optimal policy chooses the corner point &o- Therefore, 

' ^)+Q a (| 

g 2 i " * 2 



■ : , Qi(^)+Q 2 (l-D >1 e (1-6) 



e 2 

= 1 - — > 0.9700. 
Discrepant Region 2: (1+ ^~f ] < §7 < fff 

In this case the Myopic policy chooses the corner point b 2 whereas the optimal policy chooses the corner point b\. Therefore, 



8(2-e) 1 QiVs _£(£zfl_ > 9002 



4 ^ Qi V2 
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This is a minimization of a function of two variables for all possible e values in the interval < e < e t , and the ratio ^ in 
the interval (1 + e ~ e2) < < 

l-£ ^ e 

Discrepant Region 1: < ^ < 

In this case the Myopic policy chooses the corner point b\ whereas the optimal policy chooses the corner point 6q. Therefore, 



,t,/ Qii^l+MtJl > i £ , (l^il! e 

Q 2 ± 



: 2 

= 1 - — > 0.9500 



Discrepant Region 2: ^ £ e ^ < q7 < f^f 

In this case the Myopic policy chooses the corner point 62 whereas the optimal policy chooses the corner point 60 ■ Therefore, 

Qld _ I + 8(2-e) ) _ 8(2-e) ) 



>(}—±)(l- e+ e - ) + -- -! > 0.9150. 
- V 2-e A 4 4(2-e) ; 4 4(2 - e) ~ 

Discrepant Region 3: {1+ ^zf ] < §7 < 

In this case the Myopic policy chooses the corner point 62 whereas the optimal policy chooses the corner point b\. Therefore, 

Qlin — % + WTZ7\) + Q2( 3 



^, _ - s8 2 t 8 (2-e)/ ^ ^V8 8(2-e) / 

Q 1 (^)+Q 2 (i-f) 

3 e 1 e 1 Q2 ( 3 



> - 2 * ( ^ 1 8(2 - £W > 0-9474. 

4 Qi V2 4/ 

A. Weighted Departure-Rate Ratio Analysis, Case 2: e c < e < 0.5 

Considering the mappings in figures [8] and [6] again for the regions where the Myopic policy and the optimal policy "chooses" 
the same corner point, we have \P' = 1. We analyze the ratio in the regions where the two policies chooses different corner 
points termed as "discrepant" cases. Note that (1 — e)(3 — 2e) is always less than or equal to (1 — e)/e for e > e c . Since due 
to e > e c we also have < jEr, there is only one discrepancy region. 
Discrepant Region 1: (1 - e)(3 - 2e) < ^2 < 1=1 

In this case the Myopic policy chooses the corner point b\ whereas the optimal policy chooses the corner point bo. Therefore, 



^(l 2 + 8(2-e)) + Q 2 (l 8(2-e), 



> (_!_)(^_ e+ 6 - ) + -- / 6 > 0.914. 
- V l-e A 4 4(2-e) y 4 4(2 - e) - 



Combining all the cases, for all e £ [0, 0.5], we have that \L r ' > 0.9002 for all possible Q\ and Q2. ■ 
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